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1. Introduction — The Present State of the Universe 

The contemporary classical cosmology successfully describes the main 
features and evolution of the universe, but uses for this some specific initial 
conditions. In the framework of classical cosmology these conditions do not 
have their own reasonable explanation. They are just selected in such a way 
that the theoretical predictions be compatible with the actual observations. 
A more deeper understanding of why the universe has these and not the other 
properties can be provided by quantum cosmology. The most important 
unsolved issue is the nature of the cosmological singularity whose existence 
follows from the classical general relativity. The phenomenon of singularity 
is probably the most compelling reason for replacing classical cosmology with 
quantum cosmology. 

Let us recall some properties of the present world which seem to have 
their origin in the very early universe (see standard textbooks [1]). 

The distribution of galaxies in space as well as the distribution of their 
red shifts indicate that at the largest accessible scales the universe is more 
or less homogeneous and isotropic. The most convincing manifestation of 
the large scale homogeneity and isotropy of the universe is the very low level 
of the angular variations of the temperature of the microwave background 
radiation: AT/T^5 • 1CT 6 [2]. All the observational data point out to the 
conclusion that the overall structure and dynamics of the observable part of 
the universe can be described, up to small perturbations, by the Friedmann 
(or Friedmann-Robertson- Walker, FRW) metrics: 

ds 2 = c 2 dt 2 - a 2 (t)dl 2 . (1) 

It is known that dl 2 , the spatial part of the metrics (1) , can correspond 
to the closed (k = +1), open {k — — 1) or flat (k = 0) 3-dimensional spaces. 
The actual sign of the spatial curvature depends on the ratio Q = p m j p c of 
the mean total matter density p m to the critical density p c = q^qH 2 , where 
H is the Hubble parameter H(t), H(i) = a /a. The majority of the available 
astronomical data favour the value Q<1 which implies that k = — 1. However 
the observations can not presently exclude neither k = nor k — +1 . In 
any case, the current value of the parameter Q is close to one. 

Although the overall structure of the universe is homogeneous and isotropic, 
it is obviously inhomogeneous and anisotropic at scales characteristic for 
galaxies and their clusters. It is believed that these inhomogeneities were 



3 



formed as a result of growth of small initial (primordial) perturbations. In 
order to produce the observed inhomogeneities the initial perturbations must 
have had the specific amplitude and specific spectrum. There are some 
theoretical and observational arguments in support of the so-called "flat" 
Harrison-Zeldovich spectrum [3] of the initial fluctuations. In order to be 
compatible with the observations this picture may also require a significant 
amount of "dark" matter. 

The dynamical characteristics of the averaged distribution of the matter, 
the growth and formation of the small scale inhomogeneities, the abundances 
of various chemical elements, as well as other features of the actual universe, 
are all successfully brought together by the "standard" classical cosmologi- 
cal theory. The trouble is, however, that the "standard" theory postulates 
certain properties of the universe rather than derives them from more funda- 
mental principles. For instance, the observational fact of the angular unifor- 
mity of the temperature of the microwave background radiation over the sky 
does not have any rational explanation except of being a consequence of the 
postulated, everlasting homogeneity and isotropy (plus small perturbations). 
A more natural explanation to a set of observational facts can be provided by 
the inflationary hypothesis [4]. Of course, this hypothesis has its own limits 
of applicability and conditions of realization, and after all it may prove to be 
wrong, but the phenomenon of inflation seems to be quite general and stable. 
This is why it is worthwhile to investigate its consequences and compare with 
observations. 

According to the inflationary hypothesis the spatial volume of the uni- 
verse confined to the current Hubble distance In = c/Hk-2 ■ 10 28 cm or, 
possibly, even much larger volume, has developed from a small region which 
was causally connected in the very distant past. If the inflationary stage 
in the evolution of the very early universe did really take place, the large 
scale homogeneity and isotropy, as well as the closeness of Vl to one, can be 
explained as the consequences of the inflationary expansion. 

The simplest model for the inflationary stage of expansion is provided by 
the De-Sitter solution. Originally it was derived as a solution to the vacuum 
Einstein equations with a constant cosmological A-term. However, it can 
also be treated as a solution to the Einstein equations with matter satisfying 
the effective equation of state p = -e. The De-Sitter solution describes a 
space-time with a constant 4-curvature. This space-time is as symmetric 
as Minkowski space-time, in the sense that it also admits the 10-parameter 
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group of motions. The line element of the De-Sitter space-time has the form 
(see, for example, [5]): 



where a(t) = r ch^, and r = const. It is known that in the De-Sitter 
space-time one can also introduce the frames of reference with flat or open 
(hyperbolic) space sections, though these coordinate systems do not cover 
the whole of the De-Sitter space-time. The often used is the flat 3-space 
representation: 



where a(t) = e Hot , and H = const. The constant Hq is the Hubble pa- 
rameter at the De-Sitter stage of expansion. The scale factor a(t) of the 
line element (2) approaches the behaviour a(t) = e H{)t very quickly during 
several characteristic time intervals t = \ r . 

In order to see the advantages of the inflationary hypothesis let us assume 
that intially, during the Planck era, the distance between two idealized phys- 
ical objects was of order of a few Planckian scales, l v \ = 10~ 33 cm. It can be 
shown that the present day distance between these objects can be as large as 
the present day Hubble distance l H ^2 • 10 28 cm, if the duration of inflation- 
ary stage At was sufficiently long, H At > 65 In this way a small causally 
connected region could have expanded to the size of the presently observed 
universe, and the uniformity of T over the sky could have been established 
by causal physical processes in the very early universe. At the same time, a 
sufficiently long inflation makes the present value of the Q parameter very 
close to Q = I. 

2. What Can We Expect From a Complete Cosmological 



The hypothesis of the inflationary stage helps us to make some of the 
cosmological data more "natural". However, the origin of the inflationary 
stage still needs to be explained. The question still remains, what kind of 
evolution did the universe experience before the inflationary stage and how 
did the universe itself originate? A frequently made assumption is that prior 
to the De-Sitter stage there was a preceding radiation-dominated era. This 
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assumption just postpones the answer to the above mentioned questions and 
returns us to the problem of cosmological singularity and quantum gravity. 
As a more fundamental solution to the problem, it was suggested [6] that 
the inflationary era was preceded by an essentially quantum-gravitational 
phenomenon called a spontaneous birth of the universe. A theory capable 
of describing the classical stages of evolution of the universe, as well as its 
quantum-gravitational origin, can be named a complete cosmological theory. 
Let us speculate on the main expected features of such a theory. 

The desired evolution of the scale factor a(t) is shown in Fig. 1. Ac- 
cording to this scenario, the moment of appearance of the classical universe 
corresponded to t — 0. After that moment of time the inflationary evolution 
has started and has been governed by Eq. (2). It is natural to expect that all 
the characteristic parameters of the newly born universe were of order of the 
Planckian scales, i.e. the classical space-time came into being near the limit 
of applicability of classical general relativity. The inflationary expansion may 
be able to pickup such a micro-universe and to increase its size up, at least, to 
the present day Hubble radius. The wiggly line joining the points a = and 
a = l p i at Fig. 1 was meant to describe an essentially quantum-gravitational 
process similar to the quantum tunnelling or quantum decay which could 
have resulted in the nucleation of the universe in the state of classical De- 
Sitter expansion. It is reasonable to suppose that at the beginning of classical 
evolution the deviations from the highly symmetric De-Sitter solution were 
negligibly small. Moreover, it seems to be sufficient to take these deviations 
with the minimally possible amplitude, i.e. at the level of quantum zero- 
point fluctuations. During the inflationary period these fluctuations could 
have been amplified and produce the density perturbations and gravitational 
waves. The density perturbations are needed to form the observed inhomo- 
geneities in the universe. Gravitational waves seem to be the only source of 
impartial information about the inflationary epoch and the quantum birth 
of the universe. 

These matters have been subjects of study in many research papers by 
various authors. We will discuss them in more detail below. We will see that 
some of the notions introduced above have acquired more precise formulation 
and some of the problems have been partially solved. 

3. An Overview of Quantum Effects in Cosmology 

From this brief exposition of a complete cosmological theory it is clear 
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that quantum effects and quantum concepts should play a decisive role in 
different contexts and at different levels of approximation. It is useful to give 
a short classification of the areas of further discussion where the quantum 
notions will be dealt with. It is worth emphasizing that we will often use 
below the common and powerful technique which is the splitting up of a given 
problem into the "background" and "perturbational" parts. 

We will start from a description of classical perturbations on a classical 
background space-time. The physical meaning of such effects as parametric 
amplification of cosmological perturbations, and first of all, amplification of 
gravitational waves, can be clearly seen already at this level of approxima- 
tion. The next level of approximation treats the perturbations as quantized 
fields interacting with the variable gravitational field of the nonstationary 
universe, or, in geometrical language, with the classical background geome- 
try. A particular, but not obligatory, example of the variable gravitational 
field is provided by the inflationary expansion. At this level of approxima- 
tion, we will see how the initial vacuum state of the quantized fields evolves 
into a pure multiparticle state with very specific quantum properties. It will 
be shown that the final quantum state belongs to the class of the so-called 
squeezed quantum states. Squeezing is a very distinct feature potentially al- 
lowing to prove or disprove the quantum origin of the primordial cosmological 
perturbations. 

At a still deeper level, the background geometry and matter fields are also 
treated quantum-mechanically — this is the realm of quantum cosmology. 
The main object of interest in quantum cosmology is the wave function of 
the universe which, in general, describes all degrees of freedom at the equal 
footing. This level of discussion is appropriate for tackling such issues as 
the beginning and the end of classical evolution as well as quantum birth of 
the universe. However, there is no one unique wave function of the universe, 
there are many of them. All possible wave functions constitute the whole 
space of the wave functions. Presently, we do not know any guiding principle 
allowing to prefer one cosmological wave function over others. This is why 
we are facing a painful job of analyzing all of them trying to introduce a 
probability measure in the space of all wave functions. 

Going still further, one can introduce a notion of a Wave Function given 
in the space of all possible wave functions. In other words, a wave function 
of the universe becomes an operator acting on the Wave Function describing 
the many universes system. This is a subject of the now popular so-called 
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third-quantized theory. It is aimed at describing the multiple production 
and annihilation of the baby-universes. This fascinating subject is still at 
the beginning of its development and is beyond the scope of the present 
papers. The reader is referred to the recent review and technical papers on 
the subject [7]. In some sense the different theories listed from above to 
the bottom are various approximation to the theories listed in the opposite 
direction. 

Let us start from the classical theory of small perturbations superimposed 
on a given background space-time. 

4. Parametric (Superadiabatic) Amplification of ClassicalWaves 

Let us consider classical weak gravitational waves. The main purpose is 
to study the parametric (superadiabatic) amplification of gravitational waves 
[8, 9]. The same mechanism is applicable to other fluctuations if they are 
governed by similar equations. 

We assume that the space-time metric g Mi , can be presented in the form 
d^y^df} + fyui/j where gffi is the background metric: 

ds 2 = a 2 ( V ) (drf - dx 2 - dy 2 - dz 2 ) (4) 

and are gravitational wave perturbations. The functions /^(x, 77) can 
be simplified and some of them put to zero, h ^ = 0, by using the available 
gauge freedom. The remaining components can be decomposed into the mode 
functions, so that for a given mode one has 

h i = -Kv)Gi(x,y,z), (5) 

(X 

where h\ = h^g^^. 

In the case under discussion the eigenfunctions G\ can be taken in the 
simplest form: G\ = exp [ ±i(riix + n 2 y + n 3 z)} , n\ + n 2 , + n\ = n 2 . The main 
equation to be solved is 

^" + fi(n 2 -U(r ] ))=0, (6) 

where ' = d/drj, U(rj) = a" /a, n is the wave-number, and the wavelength is 
A = 2ira/n. 

Equation (6) describes an oscillator with the varying frequency, that is 
we are dealing with a parametrically excited oscillator. It is also worth 
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noting that this equation is mathematically similar to the spatial part of the 
Schrodinger equation, with the prime playing the role of a spatial derivative, 
for a particle with energy n 2 and the potential U{rj). A typical potential 
barrier U(i]) is shown in Fig. 2. 

In the intervals of ?y-time such that n 2 >>\U(r))\ the solutions to Eq. (6) 
have the form n = e ±mT? , so that one has usual high-frequency waves with the 
adiabatically changing amplitude: h = ^ sin (nrj + 0) . In the expanding 
world the amplitude decreases. The amplitudes of the waves with such n 
that n 2 »\U(r})\ for all values of i] decrease adiabatically for all rj. These 
waves are shown symbolically by the wavy line with the decreasing amplitude 
above the potential barrier U (rj) in Fig. 2. 

If for a given n there is an interval of time when n 2 <<\U(i])\, the solu- 
tions to the second-order differential equation (6) are no longer oscillatory. In 
case U{rj) = a /a they are /ii — a and /i 2 — a> J a~ 2 drj. The waves satisfying 
n 2 <<\U(r])\ for some 77 encounter the potential barrier and are governed by 
the solutions fii and fi 2 in the under-barrier region. The amplitude /i/ of the 
function f/,(r)) right after exit of the wave from under the barrier depends on 
the initial phase of the wave 0. The exiting amplitude /if can be larger or 
smaller than the entering amplitude /Xj defined right before the wave encoun- 
tered the barrier. However, averaging of (/U/) 2 over the initial phase <p (i.e., 
integrating from to 2tt ) leads always to the dominant contribution from 
the solution /ii. This means that the adiabatic factor 1 / a is cancelled out 
by /ii — a and the amplitude h (with the factor 1/a taken into account) of a 
"typical" wave can be regarded as remaining constant in the region occupied 
by the barrier. It stays constant instead of diminishing adiabatically as the 
waves above the potential barrier do. Thus, the exiting amplitude hf of a 
"typical" wave is equal to the entering amplitude hi and is larger than it 
would have been, if the wave behaved adiabatically (see Fig. 2). 

The amplification coefficient R(n) for a given n is just the ratio a(rjf) / a(^) 
where a(r)i) is the value of the scale factor at the last oscillation of the wave 
before entering the under-barrier region, and a(rjf) is the value of the scale 
factor at the first oscillation of the wave after leaving the under-barrier re- 
gion. It is seen from Fig. 2 that different waves, that is waves with different 
wave numbers n, stay under the potential barrier for different intervals of 
time. This means that, in general, the amplification coefficient depends on 
n: R(n) = 1 for all n above the top of the potential, and R(n)^>l for 
smaller n. The initial spectrum of the waves h(n) = A(n)/a, defined at 
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some 7] well before the interaction began, transforms into the final spectrum 
h(n) = B(n)/a, defined at some 77 well after the interaction completed. The 
transformation occurs according to the rule: B(n) = R(n)A(n). This is the 
essence of the mechanism of the superadiabatic (parametric) amplification 
of gravitational waves and, in fact, of any other fluctuations obeying similar 
equations. 

The initial amplitudes and spectrum of classical waves can be arbitrary. 
It is only important to have a nonzero initial amplitude, otherwise the fi- 
nal amplitude will also be zero. Now, remaining at the same classical level, 
we will imitate the quantum zero-point fluctuations by assuming that they 
are classical waves with certain amplitudes and arbitrary phases. Waves 
with different frequencies have different amplitudes, so they form some ini- 
tial vacuum spectrum. In order to derive the vacuum spectrum, we neglect 
the interaction with the gravitational field and consider, essentially, waves 
in Minkowski spacetime. The energy density of gravitational waves scales 
as e g = (c A /G)(h 2 / A 2 ). For a given wavelength A we want to have "a half 
of the quantum" in each mode, that is we want to have energy \ Hcu in the 
volume = A 3 . It follows from this requirement that the vacuum amplitude of 
gravitational waves with the wavelength A is equal to h(\) = l p i/\. Hence, 
the initial vacuum spectrum of gravitational waves, defined at some early 
epoch 77b, is /i(n)~n/ p ;/afe, where is the scale factor at that epoch. This is 
the spectrum to be transformed by the interaction with the external gravita- 
tional field. The amplification process makes the number of quanta in each 
interacting mode much larger than 1/2. The renormalization (subtraction of 
^hu ) does not practically change the energy of the amplified waves, while it 
cancels the initial \Tiuj of those (high-frequency) waves that did not interact 
with the field. 

5. Graviton Creation in the Inflationary Universe 

As an illustration, we will apply the above considerations to the inflation- 
ary (De-Sitter) model. In terms of 77 -time the De-Sitter solution (3) has the 
scale factor a{rj) = —c/H r]. (It is convenient to have 77 negative and growing 
from —00.) We assume that the De-Sitter stage ends at some 77 = 77i<0 and 
goes over into the radiation-dominated stage with 0(77) = (77 — 2r]i)c/H r]f. 
The relevant potential U{rf) is shown in Fig. 3 by a solid line 1. A wave 
with the wave number n, (nt]i) 2 <^il, enters the potential and ceases to os- 
cillate at some 77^ when A~cd/a, that is 27ra/n~a 2 /a / . This leads to the 
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entering condition nrji~l. For different n 's this condition is satisfied at dif- 
ferent a(r)i) : a(r)i)~cn/H . The waves leave the potential and start oscillat- 
ing again at the radiation-dominated stage when the condition 27ra/n~a 2 /a' 
is satisfied again. This leads to the exiting condition n(i]f — 2r]i)~l and 
a(r)f)^c/nH r)l^a(r)i)(nr)i)~ 2 . The amplification coefficient R(n) = (nrji)~ 2 
is much larger than 1 for (ni]i) 2 <^l and scales as R(n) = n~ 2 . 

Now we will see how the vacuum spectrum is transformed. By the time 
of entering the barrier, the amplitude h(rji)~n(l p i/ab)(ab/a(rji))~l p iH /c is 
the same for all n. The exiting amplitude is also independent of n : h(r]f) = 
IpiHo/c, that is all amplified waves start oscillating with the same amplitude. 
The fluctuations which start oscillating (enter the Hubble radius) with the 
same amplitude are said to have the "flat" Harrison-Zeldovich spectrum. 
One should keep in mind, however, that the comparatively shorter waves 
start oscillating earlier and their amplitudes decrease more by some fixed 
(present) time rj : h(Tj )^h(nf)(a(nf)/a )^l p in/(nni) 2 ao. As expected, the 



amplification coefficient R(n) 



n 



transforms the initial vacuum spectrum 



h(n) = n into the final spectrum h(n) = n . 

To put some more detail in this discussion above, one can consider ex- 
act solutions to Eq. (6) at the inflationary (i) and radiation-dominated (e) 
stages. In the region r\<r\\ the general solution to Eq. (6) has the form 



Hi = A 



cos (nn + 0) sin (nn + (f>) 

nn 



(7) 



where A and are arbitrary constants. For rj>rji the general solution has 
the form 

fi e = B sin (nn + x), (8) 

where B and x are constants to be determined from the conditions that \i 
and // join continuously at r\ — r\\. 

It is convenient to introduce nr\\ = x. The joining conditions are: 



A 
A 



cos (x + <f>) sin (x + 0) 

— cos (x + (p) - 1 

x V x l 



B sin (x + x) 



sin (x + 0) 



= B cos (x + x) 



(9) 



Their consequence is (B/A) 2 = 1+x 2 +(x 4 — 2x 2 ) sin 2 (x + 0) —x 3 sin 2(x + 
One can see that (B/A) 2 depends on the initial phase 0. After averaging over 



11 



the phase <p one obtains (B/A) 2 = 1 + l/2x 4 . This expression can be taken as 
the definition of the amplification coefficient R(n). It is seen that R{n) = 1 
for waves above the potential, x 2 ^>l, and 

do) 

for the longer waves, x 2 <^Cl. This expression is in full agreement with the 
previous qualitative estimates. 

The radiation-dominated era ends at some time f] 2 and goes over into the 
matter-dominated era with the scale factor 

a(rj) = c(rj + i] 2 - 4r)i) 2 /4:H r]f (r) 2 - 2r)i). For r)>r) 2 the potential U(rj) is 
again non-zero, it is shown by the solid line 2 in Fig. 3. The waves which 
satisfy the condition n(i]2 — 2r]i)^>l do not interact with the potential and 
their evolution was fully described above. (This part of the spectrum was first 
discussed in Ref. 10.) However, the longer waves, satisfying the opposite con- 
dition n{r\2 — 277i)<^l, interact with the second barrier and transform addi- 
tionally their spectrum. (This part of the spectrum was first discussed in Ref. 
11, see also Ref. 12.) These waves start oscillating at the matter-domianted 
stage when the requirement n(i]f + r\ 2 — 4?7i)~l is satisfired. For these waves, 
the exiting value of the scale factor is a(r)f)xic/ H (nr]i) 2 (r] 2 — 2r]i) and the 
amplification coefficient i?(n)~l / {ni]i) 2 n(i] 2 — 2r)\) = n~ 3 . The present-day 
amplitudes are h(j] )~l p in/ {ni]i) 2 n(i] 2 — 2i]i)a = n~ 2 . 

The waves with the wave number n# satisfying nn{vo + V2 — 4?7i)~l en- 
ter the Hubble radius at the present epoch. Their wavelengths A# are 
of the order of the present Hubble radius, \h~2 ■ 10 28 cm, and their fre- 
quencies v H are z/f/^lCT 18 Hz. The amplitude of these waves is, roughly, 
huivo^lpiHo/ - It can not be too large in order not to cause too large angu- 
lar (quadrupole) variations in the temperature of the microwave background 
radiation, AT/T^hnivo)- This requirement places a limit on the value of 
the Hubble parameter H Q at the De-Sitter stage [11]. Equally strong limit 
follows from the observational restrictions on AT/T in a few degrees angu- 
lar scale, where the contribution of waves with frequencies v = 10" 16 Hz is 
dominant [9] . The inflationary spectrum of graviational waves in terms of 
the spectral flux density as a function of frequency is, is shown in Fig. 4 
(adopted from [9]). The upper position of the spectrum is determined by 
the observational data on AT/T. At the same graph one can see theoretical 
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predictions of some other models and the existing experimental limits as well 
as the expected levels of sensitivity of various observational techniques. 

The predicted gravity-wave spectrum is more complex if the Hubble pa- 
rameter at the inflationary stage was not constant. It is interesting to know 
that the time variations of the Hubble parameter are in one-to-one corre- 
spondence with the frequency variations of the present-day spectral energy 
density of waves [13]. This makes it possible, at least in principle, to study 
the details of the very early evolution of the universe by measuring the spec- 
tral properties of relic gravitational waves. 

In our previous discussion we have been mainly interested in the ampli- 
tude of a "typical" wave. In other words, we have been calculating the r.m.s. 
value of the final amplitude, assuming that the initial phase (f> is distributed 
randomly and evenly in the interval form to 2n. However, the distribu- 
tion of the final phase is also important. Equivalently, one can ask about 
the final distributions of the quadrature components of the wave, that is the 
components proportional to sinra? and cosnrj. This study can serve as an 
introduction to the notion of the quantum mechanical squeezing which we 
will be discussing later. 

Let us return to the exact solutions (7) and (8). Initially, for r/^ — 
oo, solution (7) can be written as //j~A[i>i sin nn + v 2 cos nn] , where v\ = 
— sin0, t>2 = cos0. The mean values of v\, V2 and V1V2 are zero. However, 



the mean values of v\ and v\ are nonzero and equal: 



To derive 



the quadrature components of solution (8) one can first find, from the joining 
conditions, the constant x an d rewrite Eq. (8) in the form: 



r 1 1 

jj, e — A sin nn +— cos x sin (x + <f>) cos <p — sin 

IX 2 X 

A cos nn - sin x sin (x + </>) sin <fi + cos <fi 



+ 



x" 



X 



(11) 



Aki sin nr] + A/c 2 cos nn 

The mean values of hi and k 2 are again zero, but for the quadratic combina- 
tions one obtains: 

1 



hi 



1 




1 






f 1 




~ 2 


1 + 


X 2 


+ 


rr 2 ' 


\X 2 




1 




1 










~ 2 


1 + 


X 2 


+ 


X 2 


Kx 2 


-0 



cos 2 X 



sin 2x 



X" 



sin' 



kik 2 = - 
2 



2x 2 V 



^ x 2 



x H — - sin 2x 

x 6 



(12) 



sin 2x cos 2x 

x 6 
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Note that k\ k\ - (jbife) 



v 2 



If one makes a constant shift 
i] = fj + y, the solution (11) transforms to fi e = All sin w/ + A/ 2 cos nfj. The 
constant y can be chosen in such a way that Zi/ 2 = . Under this choice one 
obtains 





(13) 



This choice minimizes one of the quadrature variances and maximizes the 
other. 

We are interested in the case x 2 <^l. For this case, l 2 ^l/2x 4 , l\^x 4 jl. 
We see that during the amplification process, one of the noise components 
strongly increased while the other decreased equally strongly. One can also 
say that the final phase x is n °t evenly distributed as the function of (f>, but 
is highly peaked near the values tg x~ ~ % x - We will see below that the su- 
pression (squeezing) of variances in one of two quadrature components of the 
wave field is a characteristic feature of squeezed quantum states. Moreover, 
squeezing may reduce one of the two variances below the level of zero-point 
quantum fluctuations. 



6. Quantum States of a Harmonic Oscillator 

Let us first recall some properties of quantum states of an ordinary har- 
monic oscillator. We will need this information in our further discussion. 
Especially, we will be interested in the notion of squeezed quantum states. 

Classical equations of motion for a harmonic oscillator, x + oj 2 x = 0, 
can be derived from the Lagrange function L = ^mx 2 — \mw 2 x 2 according 

to the rule: (jj^) — = 0. Associated with L is the Hamilton function 
H — ^ + r ^w 2 x 2 ) where p = (§§)■ Quantization is achieved by introducing 

the operators x and p = —ih-^ and establishing the commutation relation: 
[xp] = ih. From x and p one can construct the creation and annihilation oper- 

atorsa+anda: a+ = \^)'\x - a = {^)'\x + [a, a+] = 

1 and the particle number operator N: N = a + a = -r- — \. The oscillator 
can be described by the wave function (or state function) ip(x,t) which sat- 
isfy the Schrodinger equation: = Hip. The ground (vacuum) quantum 
state |0> is defined by the requirement a|0>=0. The ground wave function 

has the form ip(x) = [^|] ^ e~^ x2 . The n-quantum states are defined as 
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the eigenstates of the N operator: N\n> = n\n>, they are also eigenstates 
of H with eigenvalues hu(n + H n> = hw(n + |) n>. These states are 
produced by the action of the creation operator a + on the vacuum state: 

| n>= (^T| >. 

An important class of quantum states are coherent states which are re- 
garded as the "most classical". The coherent states are generated from the 
vacuum state |0> by the action of the displacement operator: D(a,a) = 
exp [aa + — a*a], where a is an arbitrary complex number. We are mostly 
interested in squeezed states. Squeezed states involve operators quadratic 
in a, a + . A one-mode squeezed state (for a review see, for example, [13]) is 
generated by the action of the squeeze operator: 

Si(r,(f)) = exp ^r{e~ 2t< ^a 2 — e 2t ^a +2 ) , where the real numbers r and are 
known as the squeeze factor and squeeze angle: 0<r<oo, — | <0<|. A 
squeezed state can be generated by the action of the squeeze operator on 
any coherent state and, in particular, on the vacuum state, in which case 
it is called squeezed vacuum state. If a harmonic oscillator is exposed to 
the time-dependent interaction, a one-mode squeezed state is produced, as a 
result of the Schrodinger evolution, by the interaction Hamiltonian 

H{t) = a(t)a +2 + a*(t)a 2 (14) 

where a is arbitrary function of time. 

The meaning of the word "squeezing " is related to the properties of these 
states with respect to variances (or noise moments) AA of different operators 
A: AA = A — <A>. The squeezed wave functions are always Gaussian: 

i){x) ~ e-h* 2 (15) 

but, the variances of the variables x and p are substantially different. They 
can be presented in terms of the complex parameter 7, 7 = 71 + 272, or real 
parameters r, 0: 

1 I'vl 2 

<(Axf> = — , <(Apf> 



2 7 i ' V ' ' 2 7 i ' 
< ( Ax 2 ) > = - [chlr - shir cos 20) , ( 16 ) 



2 

9 1 

<(Ap) > = -(c/i2r + s/i2rcos20) 
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These variances should be compared with those for a coherent state, in 
which case they are always equal to each other and are minimally possible: 
<(Axf> = <(Ap) 2 > = \. So, in a squeezed state, one component of the 
noise is always large but another is "squeezed" and can be smaller than \. 
In (x,p) plane the line of a total noise K — \ <(Ax) 2 > + <(Ap) 2 > for the 
coherent states can be described by a circle, while this line is an ellipse for 
the squeezed states (see Fig. 5). In a squeezed vacuum state the mean values 
of x and p are zero, so in this case, the center of the ellipse is at the origin. 

The mean value of the particle number operator N is not zero in a 
squeezed vacuum state, it can be expressed in terms of the squeeze parameter 
r : <N> = sh 2 r. For strongly squezed states, r^>l, the <N> is very large 
<N>^\e 2r . The variance of N is <(AiV) 2 > = \sh 2 2r. The <(AiV) 2 > is 

also very large for r>l : <(AA^) 2 >^|e 4r and <(AA r ) 2 >^«^e 2r , that 

is <(AN) 2 >^<N>. In contrast, <(AiV) 2 >5 = <N>^ for coherent 
states. 

A notion which is very useful for our problem of the pair particle creation 
is the two-mode squeezed states. The two modes under our consideration 
will be two particles (waves) travelling in the opposite direction. A two- 
mode squeezed vacuum state is generated from the vacuum |0,0> by the 
action of the two-mode squeeze operator 

S{r, <f)) = exp[r(V 2i Va_ - e 2 ^a+a^)] (17) 

where a+, a_ and a\, at are annihilation and creation operators for the two 
modes, r is the squeeze parameter and is the squeeze angle. 



7. Squeezed quantum states of relic gravitons and 
primordial density pertrubations 

Our preceding analyses of perturbations interacting with the variable 
gravitational field was essentially classical. Quantum mechanics entered our 
calculations only as a motivation for choosing the particular initial ampli- 
tudes and for making the averaging over initial phases. We interpreted the 
final results as quantum-mechanical generation of gravitational waves and, 
possibly, other perturbations, but a rigorous treatment is still needed. We 
will see now that a consistent quantum-mechanical theory confirms our main 
results and gives a much more detailed and informative picture of the entire 
phenomenon. 
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We will start again from gravitational waves. The gravitational wave field 
hij(r), x) becomes an operator and can be written in the general form 



(18) 



This expression requires some explanation. In Eq. (18) we do not write the 
scale factor a{rf) in front of the expression (compare with Eq. (5)) which can 
be taken care of later. It is precisely hij 's given by Eq. (18) that appear 
automatically in the "field-theoretical" treatment of the problem, see Refs. 
[15, 16]. The normalization constant C includes all the numerical coefficients 
but we do not need them now and will not write C below. We will also 
use units c = 1, H — 1. Two tensors p^(n), s = 1, 2, represent two inde- 
pendent polarization states of each mode (wave). The tensors p\j satisfy the 
"transverse-traceless" conditions: p\p? = 0, p\^ = 0. For a wave travel- 
ling in the firection ^ = (sin 6* cosy?, sin # sin y?, cos#) the two polarization 
tensors are p^(n) = IJj — miirij, pfj(n) = krrij + Ijirii where lj, rrij are two 
unit vectors ortogonal to - and to each other: lj = (siny?, — cosy?, 0), rrij = 
(cos cosy?, cos 9 sin tp, — sin 9). The operators a*(r/), a„ + ( 7 ?) are annihilation 
and creation operators for waves (particles) travelling in the direction of n. 
These are Heisenberg operators depending on time 77. 

The evolution of the operators a* (77), a^~(rj) is governed by the Heisenberg 
equations of motion for each mode n and for each polarization s: 



da n da+ 
-i[a n , H\, 



drj 



drj 



—1 



(19) 



where H is the Hamiltonian of the problem. To derive the Hamiltonian we 
will proceed as follows. 

The classical equations (6) are the Euler-Lagrange equations following 
from the Lagrangian 



12 22 a ' I / " 1 2 
H — n ji — 2—jiji +( — ]// 



The associated Hamiltonian is 



2 1 2 2 1 n a 

p + n n + 2— up 

a 



(20) 



(21) 
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where p is the canonically conjugated momentum: p = dL/dfi' = [/,' — (a'/a)^. 
In the quantum treatment, fi and p are operators satisfying the commutation 
relation p] = i. The associated annihilation and creation operators are 

In terms of 6, 5 + , the Hamiltonian (21) acquires the form 

H = nb + b + a(r])b +2 + <r* (??)& 2 (23) 

where the coupling function cr(rj) is cr(rj) = ia' /2a. The Hamiltonian (23) 
is precisely of the form of Eq. (14). In the Schrodinger picture, the quan- 
tum state (wave-function) transforms from the vacuum state to a one-mode 
squeezed vacuum state as a result of the Schrodinger evolution with the 
Hamiltonian (23). We will use the Heisenberg picture in which the operators 
change with time but the quantum state of the system remains fixed. We 
also wish to write the Hamiltonian in the form which would explicitly demon- 
strate the underlying physical phenomenon: particle creation in pairs with 
the oppositely directed momenta. The Hamiltonian to be used in Eq. (19) 
has the form 

H = na+a n + na+ n a_ n + 2a(?7)a^al n + 2a*(rj)a n a^ n (24) 

This Hamiltonian can be derived from Eq. (23) if one considers the sum of 
two Hamiltonians (23) for two modes (with the same frequency n) b\ and 62 
and introduces a n , a_ n according to the relations: 

bi-ib 2 h+ib 2 , . 

a n = — /=-, a_„ = — ^- (25) 

The descriptions based on the 6-operators corresponds to standing waves 
and one-mode squeezed states while the description based on the a-operators 
corresponds to travelling waves and two-mode squeezed states. 

By using Eq. (24) one can write the Heisenberg equations of motion (19) 
in the form 

.dd n .(X 1 ,dd n 1 .O /r>n\ 

%—— = na n + 1— aZ n , —1—. — = ™„ — 1— a_ n (2b) 
drj a drj a 
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for each s and each n. The solutions to these equations can be written as 

a n (rj) = u n (rj)a n (0) + v n (r])a± n (0) , . 

where o n (0), a„(0) are the initial values of the operators a n (r/), a+(r/) for 
some initial time and the complex functions u n (r)), v n (r)) satisfy the equa- 
tions 

.du n . .a' „ .cfo n .a' # 

i—— = nu n + i—v n , i-— = nv n + i—u n (28) 
dry a drj a 

where |w n | 2 — \v n \ 2 = 1 and u n (0) = 1, f n (0) = 0. Note that w n , t> n depend 
only on the absolute value of the vector n, not its direction. The oper- 
ators a n (0), a„(0) (Schrodinger operators) obey the commutation relations 
[a n (0), a+ (0)] = 5 3 (n — m) and so do the evolved operators: [a n (r/), a m (r/)] = 
<5 3 (n — m) . 

It follows from Eq. (28) that the function u n + v* n satisfies the equation 

K + <)" + - ^) K + <) = (29) 

which is precisely Eq. (6). The two complex functions u n , v n restricted 
by one constraint \u n \ 2 — \v n \ 2 = 1 can be parameterized by the three real 
functions r n (rj), <j> n (rj), 9 n (r}): 

u = e i9 chr, v = e- i{fi - 2 * ) shr (30) 

which are, corrrespondingly, squeeze parameter r, squeeze angle and rota- 
tion angle 6. For each n and s, they obey the equations 

r = — cos 20 
a i 

6' = n - -sin20t/ir (31) 
a 

= — n sin 20 cthr 

a 

which can be used for an explicit calculation of r, 0, 9 for a given scale factor 
a(r]). 

The squeeze parameters have been calculated [17] (by a different method) 
for a model which we have considered in Sec. 5. The model includes three 
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successive stages of expansion: De-Sitter, radiation-dominated and matter- 
dominated. It has been shown that the squeeze parameter r varies in the large 
interval from r~l to r^l20 over the spectrum of relic gravtational waves (see 
Fig. (4)). The value r~l applies to the shortest waves with the present day 
frequencies of order of Hz, r^lO 2 is the value of r attributed to the 

waves with frequencies z/«il0~ 16 Hz, and r^l20 corresponds to the waves 
with the Hubble frequencies z/^10~ 18 Hz. 

It is important to note that Eqs. (27) can be cast in the form 

a n (r]) = RSa n (0)S + R + , a+(r]) = RSa+(0)S + R + (32) 

where 

S(r, 0) = exp[r(e- 2 ^a n (0)a_ n (0) - e 2 ^a+(0)al n (0))] (33) 
is the unitary two-mode squeeze operator and 

R(0) = exp[-^(a+(0)a n (0) + al n (0)a_ n (0))] (34) 

is the unitary rotation operator. Eq. (32) demonstartes explicitly the in- 
evitable appearance of squeezing in this kind of problem. 

We will assume that the quantum state of the field is the vacuum state 
defined by the requirement a n (0)|0> = for each n and for both s. The 
values of a n (ry), a+irf) determine all the statistical properties of the field in 
the later times. The mean values of a n (rj), a+(rj) are zero: <0\a n (r])\0> = 
0, <0|a+(?7)|0> = 0. The mean values of the quadratic combinations of 
a n(v)j a ^( r ?) are n °t zero: 

<0|a n (r?)a m (77)|0> = u n (r])v m (r))5 3 (n + m) 
<0|a+(7?)a+ (r))\0> = <(?7X(r7)5 3 (n + m) . . 

<0\a n (r])a^(r])\0> = u n (i 1 )u* m (r ] )5 3 (n-m) y ' 

<0|an(^)a m (r])|0> = v^(ri)v m (r])5 3 (n - m) 

These relationships (first two) show explicitly that the waves (modes) 
with the opposite momenta are not independent but, on the contrary, are 
strongly correlated. This means that the generated field is a combination 
of standing waves [17]. Let us see how this is reflected in the correlation 
functions of the field. 
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For purposes of illustration, we will first ignore the tensorial indices in 
Eq. (18) and will consider a scalar field 

/oo . . 

d 3 n(a n (ri)e in * + at(v)e- inx ) (36) 

It is obvious that the mean value of the field is zero, <0\h(rj, x)|0> = 0, in 
every spatial point and for every moment of time. The mean value of the 
square of the field h(rj, x) is not zero and can be calculated with the help of 
Eq. (35): 



<0\h(r), x), h(r], x)|0> = Air J n 2 dn[\u n \ 2 + \v n \ 2 + u n v n + (37) 
In term of the squeeze parameters this expression can be written as 

rco 

<0\h(r], x), h(r), x)|0> = An n 2 dn(ch2r n + sh2r n ■ cos20 n ) (38) 

J o 

(it includes the vacuum term 4-7T Jq 00 n 2 dn, which should be subtracted at the 
end). It is seen from Eq. (38) that the variance of the field does not depend on 
the spatial coordinate x. The function under the integral in Eq. (38) is usually 
called the power spectrum of the field: P(n) = n 2 (ch2r n + sh2r n cos 2(f) n ). 
The important property of squeezing is that the P{n) is not a smooth func- 
tion of n but is modulated and contains many zeros or, strictly speaking, very 
deep minima. To see this, one can return to Eq. (31). For late times, that is, 
well after the completion of the amplification process, the function a' /a on 
the right-hand side of Eq. (31) can be neglected. This is equivalent to saying 
that one is considering waves that are well inside the present day Hubble ra- 
dius. For these late times, the squeeze parameter r n is not growing any more 
and the squeeze angle <fi n is just <p n = —nr\ — §§ n . Since r n ^>l for the frequen- 
cies of our interest, the P{n) can be written as P(n)~n 2 e 2r ™ cos 2 {nrj + 4>o n ). 
The factor cos 2 (nrj + On ) vanishes for a series of values of n. At these fre- 
quencies, the function P(n) goes to zero. The position of zeros on the n axis 
varies with time. Similar conclusions hold for the spatial correlation function 
<0\n(r], x), h(r), x + l)|0> : 



00 2 sin nl 



<0\h(n, x), h(n, x + 1)|0> = Air / n — — (ch2r n + sh2r n cos2(j) n )dn 

J Tbb 



(39) 
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The resulting expression (39) depends on the distance between the spatial 
points but not on their coordinates. The power spectrum of this correlation 
function is also modulated by the same factor cos 2 (nrj + 4>Q n ). It is necessary 
to note that the power spectrum of the energy density of the field is smoooth 
as it includes, in addition to Eq. (39), the kinetic energy term 

1 f°° 

-=<Q\ti(n, xWfa, x)|0> = 4vr / n 2 dn(ch2r n - sh2r n cos 20 n ) (40) 
n 2 Jo 

so the contributions with the oscillating factors cos20 n cancel out. 

Let us mention some new features which arise when one considers the 
tensor field hij(r], x), Eq. (18), not the scalar field h(r], x). First, one sees 
from Eq. (28) that the solutions u s n {rj), v'^(i]) for both polarisations s = 1, 2 
are identical; they obey the same equations with the same initial conditions. 
This means that, for each n, both polarisation states are necessarily gener- 
ated, and with equal amplitudes. This feature can serve as a clear distinction 
between gravitational waves generated quantum-mechanically and by other 
mechanisms. The second feature is related to the properties of the corre- 
lation function analogous to Eq. (39). In case of the tensor field hij, there 
is one combination of the components which has a particular interest: 
h(e) = h i je l e : >, where e l is an arbitrary unit vector. The h(e) enters the cal- 
culations of the fluctuations of the microwave background temperature seen 
in the direction e % (Sachs- Wolfe effect): 

AT / -\ 1 rm dh, H 



rM-sjC^ (41) 

The relevant correlation function is <0\h i je l e : '(r], 0)h i je t e J ' (rj, re fc )|0>, 
where r is a parameter along the line of sight. It is interesting to calcu- 
late this function for the initial time rj = and, also, for some very late time 
77. Without going into the details (they will be published elsewhere), we will 
present some results. For 77 = one has v n (r)) = 0, u n (i]) = 1 and 



g—in,je J T 



<0|/i lje V(0, 0)/^eV(0, re fc )|0> = j°° d 3 n (pl/e j f + (pJeV) 

(42) 

The presence of the both polarisations is very important: because of this 
the integration over angular variables eliminates dependence on the direction 
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e> . The final result is 

<0|/^eV(0, 0)/iyeV'(o, re fe )|0> = 16^^ J°° (nr)" 5/2 J 5/2 {nT)n 2 dn 

(43) 

One can see that the (vacuum) correlation function depends only on the 
distance between the points and does not depend on the direction e- 5 . 

To conclude this section, we should say that the same theory of squeezed 
states is applicable to the density perturbations generated quantum-mechanically 
[18]. In a simple inflationary model governed by a massive scalar field, the 
progenitor of the density perturbations and, later, the density perturbations 
themselves, satisfy equations similar to the equations for gravitational waves: 
v" + (n 2 — V(n, r)))v — where v is a gauge-invariant function which includes 
perturbations of the matter variables and gravitational field [19], and V(n, rj) 
is the effective potential analogous to the gravitational- wave potential U(rj), 
Fig. 3. Similarly to the case of gravitational waves, squeezing in density per- 
turbations and associated (longitudinal) gravitational field, exhibits distinct 
observational features. 

8. Quantum Cosmology, Minisuperspace Models and Inflation 

Until now we have been discussing the quantum fluctuations superim- 
posed on a given classical background spacetime. The next level of complex- 
ity is the quantization of the background geometry itself. This is the domain 
of full quantum gravity and quantum cosmology. 

In canonical quantum gravity the role of a generalized coordinate is played 
by a 3-geometry g( 3 \ The full set of all 3-geometries forms a superspace, 
where the wave function of the quantized gravitational field is defined. If 
some matter fields are present, the superspace includes the matter variables 
as well. In cosmological applications, one usually considers topologically 
compact geometries and calls the wave function of the entire system the 
wave function of the universe. The basic equation which governs the wave 
function of the universe is called the Wheeler - DeWitt (WD) equation. (For 
reviews of quantum gravity and quantum cosmology, see, for example [20].) 

A simplified problem, which allows a detailed investigation, is provided 
by minisuperspace models. In minisuperspace models one neglects all de- 
grees of freedom except of a few. A quantum cosmological model describing 
a homogeneous isotropic universe filled with a massive scalar field gives a rea- 
sonably simple, though sufficiently representative case. In this case one has 
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only two degrees of freedom (two minisuperspace variables): the scale factor 
a(t) and the scalar field 0(t). Since the formulation of a quantum problem 
includes integration of some quantities, such as the Hamiltonian function, 
over 3 - volume, one normally considers closed 3-sphere geometries, k = +1, 
or torus-like geometries, k=0, in order to avoid infinities arising because of 
spatial integration. The total energy of a closed world is zero. This is why 
the analog of the Schrodinger equation takes the form H-0 = 0, which is the 
Wheeler - DeWitt equation. 

For a FRW universe filled with a scalar field 0, V(<p) = ^m 2 <f) 2 , the 
Wheeler - DeWitt equation can be written as follows [21]: 

! °~ ka 2 + m 2 2 a 4 )i/j(a, 0) = (44) 




a 2 d(f) 2 



The factor p reflects some ambiguity in the choice of operator ordering. 
The possible preferred choice of p for the given model is p = 1. 

First, we will show how classical Einstein equations of motion follow from 
the quantum equation (44) in the quasi-classical approximation. For sim- 
plicity we consider the limit where the spatial curvature term ka 2 can be 
neglected. In this limit (and for p — 1) Eq. (44) reduces to 
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In the quasi-classical approximation, the wave function ip( a, 0) has the 
form ip(a,<f>) = exp(iS(a, 0) + ia(a, 0) + ...). By using this representation, 
the following equations can be derived from Eq (45): 

. d 2 S _ 2 dS_da + ?_dS__ ±^S_ + ^dS_da _ Q 
da 2 da da a da a 2 <90 2 a 2 (90 (90 

Eq. (46) is the Hamiltonian- Jacobi equation for the action S. A real 
solution to Eq. (46), which describes the classical dynamics of the model, 
can be presented in the form 

S(a,0) = -a 3 /(0) (48) 
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where the function f((f>) , as a consequence of Eq. (46), satisfies the ordinary 
differential equation 

9/2 -(^) 2 = mV (49) 

The classical equations of motion can be obtained from Eq. (48) and the 
Lagrangian L = |(— ad 2 + a 3 — a 3 m 2 2 ) of the system in the usual way. 

One writes || = —ad — || , || = a 3 = || which leads to the relations 

- = 3/, = -/' (50) 
a 

where the prime denotes a derivative with respect to 0, and the dot denotes 
a derivative with respect to the time t. By taking time derivatives of Eq. 
(50) and using Eq. (49), one derives the equations of motion which can be 
cast in the usual classical form: 



+ 3^0 + m 2 = O, f^ 2 = 2 + m 2 2 
a \aj 



(51) 



Eqs. (51) are invariant under the transformation t— > — t. The three 
equations of motion (51) can be combined in one equation, in which the time 
parameter t does not appear at all: 



da z \ da 



1- ^ 




= (52) 



(For convenience, we use the variable a = In a here and below.) This 
equation completely describes the classical trajectories in the (a, 0) space. 
The direction of motion along the trajectories is determined by the choice of 
the time direction. 

At this point we should comment on whether or not the sign of the ac- 
tion S has anything to do with the expansion or contraction of a cosmological 
model. This issue is often discussed in the context of the so-called "tunnel- 
ing" wave function [22]. Eqs. (48) and (50) may lead to the impression that 
S<0(f>0) corresponds to the expansion (d>0), while the opposite choice 
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S>0 (/<0) corresponds to the contraction (d<0) of the cosmological vol- 
ume. However, the choice of the parameter t in these equations is absolutely 
arbitrary. The functions S>0 (/<0) can perfectly well describe expansion if 
one changes the parameter t to —t in Eq. (50). Thus, the sign of the action 
does not prescribe a particular meaning to the direction of evolution along 
the classical trajectories. 

For the model defined by classical equations of motion (51), all trajectories 
of the model in the <j)j phase plane have previously been found [23]. It 
has been shown that in the case of expansion (i.e., a > 0), all the trajectories, 
except for two, start out from the ejecting nodes K 1 and K 2 (see Fig. 6). The 
remaining two trajectories, corresponsing to the inflationary regime, are two 
attracting separatrices that originate at the saddle points Si and S 2 - The 
solutions to Eq. (49) have the following asymptotic behavior for trajectories 
that start out from the nodes: f^ce ±3(t) , c 2 e ±3< ^>m 2 2 , c = const. And for 
the separatrices one has /~±|m</>, 90 2 ^>1. Different values of the constant 
c select different trajectories leaving the nodes. Hence, a particular solution 
to Eq. (49) gives a difinite function S and, at the same time, a particular 
classical trajectory. 

Now we will relate different wave functions to different classical solutions. 
One can distinguish different solutions to Eq. (49) by the subscript n which 
varies continuously and takes on two distinct values corresponding to the 
separatrices. By virtue of the linearity of the WD equation, we can present 
a full set of solutions to Eq. (45) in the form ip = J2n ex P(.iA n + iS n ), S n = 
—exp(3a)f n , A n = const valid in the lowest approximation. To every quasi- 
classical wavefunction ip n = exp(iS n ) one can put into correspondence a fam- 
ily of lines that are ortogonal to the surfaces S n = const (Fig. 7). These 
surfaces are constructed in the minisuperspace (a, 0) endowed with the met- 
ric tensor 

G* v — e~ 3a diag(-l, +1), //,z/ = l,2, x 1 = a, x 2 = <p (53) 

The vector ortogonal to S n = const can be obtained by acting on ip n 
= exp(iS n ) with the momentum operators n a and 7if 




i da 
ld_ 

i d(f) 



n 




d<f) 



n 



1p n = N^, 



n 
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Taking into account Eqs. (48), (53) one obtains N a = 3f, = . The 
vector field (iV Q , determines the lines x^(a, 0) ortogonal to S n = const 
in a parametric way: dx^/dt = N^. These lines coincide with the classical 
trajectories (see Eq. (50)). One can also note that by integrating the re- 
lation da/dcj) = N a /N^ = —3f/f ; along every classical path in the (a, 0) 
plane one gets z(a, 0) = const where z = a + 3 / (/ / f')d(f>. 

In the case at hand, the family of lines ortogonal to S n and the associated 
tangent vectors N a , are independent of a and transform into themselves 
under the shift a — * a+ const, or a(t) ^const a(t). This symmetry is a 
reflection of the fact that the function a(t) alone does not appear in Eq. (51), 
it appears only as the Hubble factor a/ a. Therefore, the invariance of the 
vector field under the displacement a — > a + const means that, for a given 
S n , the lines traced out by are all copies of one and the same physically 
distinct classical solution. It happened as a consequence of our assumption 
of a negligibly small spatial curvature, k=0; in general, it is not the case. 

Thus, we see that different solutions to the Hamilton- Jacobi equation 
determine different wave functions in their lowest (in terms of h ) approxi- 
mation. On the other hand, to a given S n one can assign a family of classical 
trajectories. The next approximation to S n defines the prefactor to the wave 
function ip n = e . The prefactor is responsible for forming a packet from 
classical trajectories determined by S n . It assigns different "weights" to dif- 
ferent classical paths. 

Let us return to Eq. (47) for a n . The general solution for a n can be 
expressed in terms of the function f n (0): a n (a, 0) = f (3a + In f'^j + B n (z) 
where B n is an arbitrary function of its argument z. In the considered ap- 
proximation, the general solution to WD equation can be written in the 
form 

^ = Yl Z ^n = Z n exp(iS n + ia n ) (54) 

n n 

where 

^ = Xne^(a"f' n y 1/2 e-^, (55) 

Xn and 7„ are arbitrary real functions of z and Z n are arbitrary complex 
numbers. One can see that to every path z = const in (a, 0) plane one can 
assign a number Q n = Xn 2 (z) which is conserved along this path. A partic- 
ular value of Q n is determined by a chosen wave function (in other words, by 
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the chosen boundary conditions for the wave function) and, specifically by 
the function Xn( z )- D iff rent wave functions favor the inflationary trajectories 
to a different degree (see, for example, [24]). 

9. From the Space of Classical Solutions to the Space of 
Wave Functions 

From the problem of distributing "weights" among different classical tra- 
jectories belonging to the same family determined by S n we now turn to the 
more difficult problem of distributing "weights" among the wavefunctions 
themselves. As we saw above, the WKB components ip n , Eq. (55), partici- 
pate in the general solution Eq. (54) with arbitrary complex coefficients Z n . 
They determine one or other choice of possible wavefunctions. How can one 
classify the space of all possible wavefunctions? 

To answer this question we will start from the simplest situation, when 
the number of the linearly independent solutions to WD equation is only 
two. For this aim we will first consider Eq. (44) in another limiting case, 
namely, when the term — -§^2 can be neglected. In this case the variable 
plays the role of a parameter and the problem reduces to a one-dimensional 
problem. The basic Eq. (44) can be written down in the form (for k = +1): 



where H = m 2 (j) 2 . We prefer to work with exact solutions to Eq. (56) so we 
choose p = -lorp = 3 [25]. (The case p = - 1 was first considered in Ref. 
[26].) We will write the exact solution for the p = - 1 case in the form: 



where u = H 2 a 2 — 1 and I1/3, K4/3, H^\/3, lT 2 \/3 are the Infeld, Macdonald 
and Hankel special functions, correspondingly. 




(56) 




(57) 
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Eq. (56) has the form of the Schrodinger equation for a 1-dimensional 
problem with the potential V(a) = a 2 - H 2 a 4 (see Fig. 8). The coefficients 
Ai, A 2 and B 2 are two pairs of arbitrary constant coefficients in front 
of two pairs of linearly independent solutions. By matching the solutions 
(57) and (58) at the point a = 1/H one finds [25]: B 1 = -Ax(l + iVz) - 

A 2 (l-iV3), B 2 = %(A 2 -A 1 ). 

Now let us characterize the full space of the wave functions. (Here we 
mainly follow Ref. [27].) In the present case, our quantun system has only 
two linearly independent states and therefore resembles a simple spin - 1/2 
system. Let us call an arbitrarily chosen basis of states |1 > and |2 >. A 
general state \ip > can be expanded as \ip> = Z x | 1> + Z 2 | 2 >, where Z 1 
and Z 2 are arbitrary complex constants. 

It is a general principle of quantum mechanics that state vectors which 
differ only by an overall non-zero multiple A describe one and the same phys- 
ical state. Thus, the pair of coordinates (Z 1 and Z 2 ) and the pair ( XZ 1 and 
AZ 2 ) are equivalent. It follows that physical quantities can only depend on 
the ratio ( = Zi/Z 2 which is invariant under rescaling. In our example above 
we may identify Z\ with Bi and Z 2 with B 2 . It is convenient to introduce 
the notation B\ = \Bi\exp(if3i), B 2 = \B 2 \exp(i(3 2 ), (5 — fli — (5 2 and then 
( = B\jB 2 = xexp{if3). The ratio ( parameterizes the points on a 2 - di- 
mensional sphere and so we see that the set of possible wavefunctions is in 
1-1 correspondence with the points on the 2 - sphere. 

We now wish to place a measure on the space of quantum states. Of 
course there are many possible measures. However, in choosing a measure 
we should be guided by the priniciple that the measure should be independent 
of the arbitrary choice of basis states |1 > and |2 >. That is if we perform a 
unitary change of basis, which will preserve all probability amplitudes, then 
the measure should remain invariant. 

The invariance of the measure may be taken as the quantum analogue of 
the principle of general covariance in classical general relativity. In fact in the 
classical limit it corresponds to invariance under canonical transformations. 
This latter invariance was used in Ref. [28] to suggest a suitable measure on 
the set of classical solutions. 

For a 2-state system the 2-dimensional unitary transformations will act 
(provided |1 > and |2 > are normalised) on the complex 2-vector (Z x and 
Z 2 ) by multiplication by a 2 by 2 unitary matrix. Clearly the ratio ( is unaf- 
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fected by matrices which are merely multiples of the unit matrix so we may 
confine attention to special unitary matrices of determinant unity, this still 
allows minus the identity matrix so if we want just the transformations which 
change the physical states we must identify to SU(2) matrices which differ 
by multiplication by minus one. That is, the effective physical transforma- 
tions acting on the space of quantum states is the rotation group SO (3) = 
SU(2)/C 2 where C 2 is the group consisting of + 1 and - 1. In fact this acts 
on the 2-sphere in the usual way provided we identify (5 with the longitudinal 
angle and x = cotan(0/2) where 9 is the usual co-latitude. 

It is now clear that we must choose for our invariant measure on the 
space of quantum states the usual volume element on the 2-sphere. This is 
clearly invariant under rotations and up to an arbitrary constant multiple it 
is unique. That is the measure in terms of (5 and 9 is: 

dV = sin9d9dp, 0<9<n, 0</3<2tt (59) 

Of course the measure is just the Riemannian volume element with respect 
to the standard round metric on the 2-sphere. 

It should be mentioned that the well known Hartle-Hawking wavefunc- 
tion [29] is exactly the south pole (9 = ir) of the 2-sphere. This wavefunction 
is real. Another real wavefunction corresponds to the north pole of the 2- 
sphere. We call this wavefunction anti-Hartle-Hawking wavefunction. All 
other wavefunctions are complex. 



10. On the Probability of Quantum Tunneling from "Nothing" 

The measure introduced in the space of all wave-functions may allow us 
to formulate and solve some physically meaningful problems. We will try to 
pose one such a problem already in the considered simplest model. As was 
mentioned above, Eq. (56) looks like the Schrodinger equation for a particle 
of a zero energy moving in the presence of the potential V(a). The form 
of the potential (Fig. 8) motivates the expectation that some of the wave- 
functions may be capable of describing the quantum tunneling or decay. In 

ordinary quantum mechanics the quantity D, where D = ^4 , can be 
interpreted as the quasiclassical probability for the particle to tunnel from 
one classically allowed region to another (see Fig. 8). The wave function 
used in this expression is determined by the imposed boundary conditions, 
i.e. it is determined by the physical formulation of the problem. The value 
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of D is always (much) less than unity, D<1, for wave functions describing 
quantum tunneling or decay. One can define a similar quantity D in our 
quantum cosmological model, though, the physical interpretation of D is 
less clear. The main difference is that in ordinary quantum mechanics one 
imposes suitable boundary conditions in time t and space x, while in our 
problem there is only one coordinate, a. (The notion of the break of classical 
evolution in quantum cosmology is rather delicate. We have argued in Ref. 
[30], that only in superspaces of more than one dimension, this notion can 
be clearly formulated.) Nevertheless, we will adopt the same definition of D 



in our problem: D = 
defined mathematical 



, , , 2 



, and will consider p<l. The quantity D is well 



V(o) 

y and can be calculated for every solution of Eq. (56) 
regardless of its interpretation. Since in quantum-mechanical problem the 
"energy" e of the particle is e = 0, we can provisionally interpret D as the 
probability for creation of the universe from "nothing" . (This is more precise 
formulation of the notion, introduced at the beginning of these lectures and 
graphically depicted in Fig. 1.) Therefore, we are interested in dividing 
the wave functions into two classes which predict D<1 and D>1. It is not 
excluded that the wave-functions predicting D<1 can be eventually justified 
as describing the quantum tunneling from "nothing", or rather from the 
"vacuum" defined in the framework of some more deep quantum theory. 

It is easy to calculate D in the approximation H << 1 [25]. One can see 
that the different choices of the wave- function give different values of D. For 
instance, the Hartle-Hawking wave-function corresponds to the choice B^O 
and gives D = expi^^^l. We are interested to know the value of D for 
a typical wave-function. In other words, we need to know how many wave- 
functions give D<1 or D>1? To answer this question one must consider the 
space of all possible wave-functions with a suitable measure. By using the 
measure (59) one can show that the set of wave functions predicting D>1 is 
very small compared with that predicting D<1. This follows from the fact 
that the surface area of the patch covered by the wave-functions with D>1 
is very small compared with the total surface area of the 2-sphere. Indeed, 
the circle separating D>1 and D<1 regions on the 2-sphere corresponds to 
the value Oq^tt — 2 exp(^— 3^2), vr — 6 <^il. Only a small area around the 
south pole 9 = n gives the wave-functions with D>1, the rest of the surface of 
the 2-sphere corresponds to the wave- functions with D<1. The ratio of the 
surface area around the south pole to the total surface area is very small; it is 
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equal to expy—j^J^l. Thus, one can say, that the probability of finding a 
wave function with D>1 (among them is the Hartle-Hawking wave-function) 
is very small. One can conclude that the overwhelming majority of the wave- 
functions seem to be capable of describing the quantum tunneling or decay, 
since they predict D<1. (It is interesting to note that the product of surface 
areas with D>1 and D<1 to their corresponding maximal values of D gives 
approximately equal numbers, both of order unity.) 

The simple example presented above clarifies the notion of the measure in 
the space of all physically distinct wave-functions. In a similar way one can 
introduce the measure in the multidimensional space of the wave-functions 
described by Eq. (54) [27]. The use of this measure shows that the inflation is 
indeed a property of a typical wave function, at least, under some additional 
assumptions adopted in [27, 25]. 

11. Duration of Inflation and Possible Remnants of the 
Preinflationary Universe 

In the framework of the inflationary hypothesis, one normally considers 
cosmological models whose period of inflation lasted much longer than the 
minimal duration necessary to increase a preinflationary scale to the size of 
the present-day Hubble radius. In such models, the number N of e-foldings 
of the scale factor during inflation is much larger than the minimal N m i n , in 
which case the volume covered by inflation is much larger than the present- 
day Hubble volume. For this reason, one normally does not expect to find any 
"remnants" of the preinflationary universe (see, however, Ref. [31]) as they 
were enormously diluted and spread over the huge inflated volume. Never- 
theless, according to the quantum cosmological considerations, the duration 
of inflation close to the minimally sufficient amount may happen to be the 
most probable prediction of some popular quantum cosmological models, as 
we will see below. (This section is based on Ref. [32]). 

Quantum cosmology is supposed to provide initial data for classical cos- 
mological models and resolve such issues as the likelihood of inflation and 
its probable duration. Obviously, we are still far away from a satisfactory 
answer. A part of the problem is that there are too many possible wave 
functions: the trouble of selecting an appropriate classical solution from the 
space of all possible classical solutions is replaced by an even bigger problem 
of selecting an appropriate wave function from the space of all possible wave 
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functions. However, if a cosmological wave function is chosen, the derivation 
of the probability distribution of the permitted classical solutions seems to 
be more straightforward. 

A wave function which has received much attention in the literature is the 
Hartle-Hawking wave function ipHH- As we have seen above, one cannot say 
that the ipHH is i n any sense more probable than others. On the contrary, 
it looks, rather, as an exception. For simple quantum cosmological models 
allowing inflation, the Hartle-Hawking wave function corresponds to a single 
point — a pole on the two-sphere representing the space of all physically 
different wave functions (see Sec. 9). However, the ipHH is a real wave function 
while all others (except the one corresponding to the opposite pole which 
is also real and which we call the "anti-Hartle-Hawking" wave function) are 
complex. This exceptional property of the i/jhh alone, if for no other reasons, 
justifies special attention to this wave function and makes it interesting to 
see what kind of predictions with regard to inflation follow from it. 

For the case of homogeneous isotropic models with the scale factor a(t) 
and a scalar field (f>(t), the ipHH predicts a set of classical inflationary so- 
lutions which can be described as trajectories in the two-dimensional space 
[a(t), 4>{t)} (see [21] and Sec. 8). These trajectories begin in the vicinity 
of a line which is the caustic line for the so-called Euclidean trajectories. 
The probability distribution Phh for the classical (Lorentzian) inflationary 
solutions follows from the ipHH and has the form 

(60) 

where N is the normalization constant and H{<p) is the Hubble factor at 
the beginning of inflation. The function Phh varies along the caustic line 
and increases rapidly toward the smaller values of <fi. This means that the 
probability to find a given inflationary solution is higher the lower the initial 
value of the scalar field (pit) (if, of course, this interpretation of Phh is 
correct). But smaller initial values of <p{t) correspond to the shorter periods 
of inflation which makes solutions with a shorter period of inflation much 
more probable than solutions with a longer period of inflation. 

An important fact is, however, that the inflationary period cannot be 
too short. The reason is that the caustic line does not extend down to the 
very low values of 0; instead, it has a sharp cusp (singularity) at the point 



Phh = Nexp 



3H 2 ((f)) 
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of return from which the second branch of the caustic line develops (see 
Fig. 9) [24]. The point of return on the caustic line divides the Euclidean 
trajectories into two families which touch the first or the second branch of 
the caustic, respectively. The Lorentzian inflationary solutions cannot begin 
with the initial value of the scalar field and the Hubble factor lower than the 
value corresponding to the point of return 0* and, therefore, their periods of 
inflation cannot be arbitrarily short. Thus, ipHH gives more weight to infla- 
tionary solutions with lower initial values of <f>(t) but does not accommodate 
solutions which begin with <f)(t) smaller than 0*. The numerical estimates 
for the case of the scalar field potentials V(0) = m 2 2 /2 and V((f>) = A0 4 /2 
show [24] that the number 0* falls short a factor 4 or 3, respectively, to ensure 
the minimally sufficient inflation. The inflated scale turns out to be of order 
10 21 cm instead of the required 10 28 cm. At the same time, the probability 
distribution function P HH reaches its maximum value at = 0*. Thus, it 
seems that the most probable prediction of the Hartle-Hawking wave func- 
tion is a "small, underinflated universe". However, it is possible that the 
discrepancy between In and the predicted inflated scale may be weakened 
or even removed for other scalar field potentials. Apart from that, the defi- 
ciency of 0* in being just a numerical factor 4 or 3 smaller than necessary, 
in the situation where the initial values of the scalar field can vary within a 
huge interval from 0* up to about 10 5 0*, can serve as an indication that the 
duration of inflation close to the minimally sufficient amount should, proba- 
bly, be taken seriously, at least, as a prediction of the Hartle-Hawking wave 
function. 

The meaning of the above discussion is that the search for the "remnants" 
of the preinflationary universe, in the framework of the inflationary hypoth- 
esis, may not necessarily be of a purely academic interest. 

12. Relic Gravitons and the Birth of the Universe 

The quantum cosmological mini-superspace models analyzed above in- 
cluded only two degrees of freedom and corresponded to homogeneous isotropic 
universes. The inclusion of all degrees of freedom at the equal footing would 
present a formidable problem. However, this problem can be simplified in 
a perturbative approximation which is a quantum-mechanical treatment of 
a perturbed homogeneous isotropic universe. In particular, the Schrodinger 
equation for gravitons, with the Hamiltonian equivalent to Eq. (21), can be 
derived from the fully quantum cosmological approach as an approximate 
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equation for the linearized perturbations. 

Let us consider a closed universe governed by an effective cosmological 
term A and perturbed by weak gravitational waves. The WD equation for 
the wave function of this system can be written in the form (see, for example, 
Ref. [33], [34]): 



2 



,3 
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(61) 

Here, the h n ; m denotes the amplitude of the gravity wave perturbation in 
a given mode (n,l,m). Since we are working in a closed 3-space, it is conve- 
nient to attribute the indices l,m to spherical harmonics. The H ra ; m denotes 
the Hamiltonian of the perturbation: H n i m = n 2 llm /2M + MO^/i^ im /2, where 
TTnim is the momentum canonically conjugated to h n i m , and M = a 3 , Q n = 
a _1 (n 2 — l) 1 ^ 2 . In what follows we will often omit the indices n,l,m for sim- 
plicity. 

The total wave function ip(a, {h}) depends on a scale factor a and a set 
of the gravity wave variables h n \ m . The ?/>(a, {h}) can be presented in the 
form ip(a, {h}) = exp[— A(a) — Ai(a)]$(a, {h}), where A(a) is the "unper- 
turbed" (background) action and Ai(a) is the prefactor of the background 
wave function. The $(a, {h}) is the part of the total wave function describ- 
ing the fluctuations. We assume that ^(a, {h}) satisfies the quasiclassical 
approximation with respect to the variable a. This allows us to simplify 
the WD-equation. We assume also that the fluctuations are weak and do not 
affect the background so that the term d 2 &/da 2 can be neglected in Eq. (61). 

It follows from Eq. (61) that the wave function ip(a,, h) for each mode of 
fluctuations obeys the Schrodinger equation 

~% = B4 (62) 

where d/adr] = —ia~ 1 (dA/da)d/da, and H = H n i m . The wave function 
<3>(a, {h}) is constructed as follows: $ = I\ n i m ^>- One can see that in the 
regime when A(a) describes classical Lorentzian evolution, that is, when the 
background space-time is the De-Sitter solution, Eq. (62) coincides with the 
Schrodinger equation for the problem considered in Sec. 7. (One has to 
take into account some obvious modifications related to the fact that we are 
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considering now the case k = +1.) However, Eq. (62) has, in fact, a wider 
domain of applicability. The assumptions under which Eq. (62) was derived 
retain this equation valid in the region where A (a) describes a classically 
forbidden behaviour of the universe, i.e. this equation is valid in the under- 
barrier region a<l/H (see Fig. 8) as well. In this region Eq. (62) takes 
the form of the Schrodinger equation written in the imaginary time. Thus, 
the graviton wave function ^(a,h) extends to the classically forbidden region 
a<l/H and may be sensitive to the form of the background wave function in 
this region. 

Our final goal is to show that the initial quantum state of gravitons at 
the beginning of the De-Sitter stage (before the parametric amplification has 
started) is not unrelated to the form of the background wave function of the 
universe in the region a<l/H. Everywhere in our previous discussion we 
were assuming that the initial state of gravitons at 77 = rjb was the vacuum. 
The present-day observational predictions have also been derived under this 
assumption. However, this assumption, though quite usual and natural, is 
not obligatory. If the initial state of gravitons could have been a non- vacuum 
state, then it would lead to the differing predictions for the present day 
spectrum of relic gravitons and their squeeze parameters. In this way, by 
measuring the actual parameters of relic gravitons, one could learn something 
about the wave function of the universe in its classically forbidden regime. 

One should note, however, that the possible deviations of the initial quan- 
tum state of gravitons from the vacuum state, regardless of the origin of 
these deviations, can not be too large. These deviations should satisfy two 
requirements. First, they should not violate our basic assumption that the 
back-action of gravitons on the background geometry is always negligibly 
small. Second, they should not lead to the predictions for the present day 
amplitudes which would exceed the existing experimental limits. By combin- 
ing these requirements one can show that only for low-frequency waves and 
only for cosmological models with minimally sufficient duration of inflation 
the initial quantum state of graviton modes can possibly deviate from the 
vacuum [35]. In this case the deviations of the present-day spectrum can be 
as large as is shown by the broken line in Fig. 10 for a specific De-Sitter 
model with l = c/ H = 10 9 l p i. In this figure the dotted line shows the 
spectrum produced from the initial vacuum state in the same model, and the 
solid line shows the highest possible inflationary spectrum compatible with 
the observational limits. (The solid line is just a low frequency part of the 
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inflationary spectrum presented in Fig. 4). 

Now we return to the question of which of the background wave functions 
are compatable with the deviations of the initial quantum state of gravitons 
from the vacuum. As we have already seen, the Hartle-Hawking wave func- 
tion ipHH and the anti-Hartle-Hawking wave function ip a HH are, in a sense, 
two extremes in the description of the classically forbidden domain. Each of 
these extremes can be used in Eq. (62) as a background wave function. For 
each of them the solution to Eq. (62) can be presented in the Gaussian form 
(compare with Eq. (15)) tp(h,rf) = C(r])e~ B ^ h . Here B(rj) is constricted as 
a linear combination of two independent complex solutions to the classical 
wave equation (6), and the vacuum state at rj = r/b corresponds to the value 

B(rj) = B(r)f,), where B{rj b ) = u b = {n 2 — l) 5 aT 2 (r] b ) . It is important 
that the choice of ipHH or tp a HH m the classically forbidden region restricts 
the function B(rj) in different ways if one is willing to subject the wave func- 
tion ip(h, rj) to the condition of normalizability: J^^ipdhKoc [34]. If this 
condition is imposed, it requires Re B(r])>0. It turns out that this require- 
ment singles out the vacuum initial state if the background wave function 
is iPhh and it leaves room for non- vacuum initial states if the background 
wave function is ip a HH- Thus, if relic gravitational waves are detected with 
properties different from those following from the initial vacuum state one 
could conclude that the universe was described by ip a HH, and not by ipHH, in 
the classically forbidden regime. This would strengthen the hypothesis that 
the universe was created in a quantum process similar to quantum tunneling 
or decay. Thus, the difference between possible wave functions of the uni- 
verse in the classically forbidden regime can be distinguished by exploring 
the properties of the gravitational wave background existing now. 
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Figure Captions 

1. The scale factor of a complete cosmological theory. 

2. Parametric (superadiabatic) amplification of waves. 

3. The potential U(t]) for the inflationary — radiation-dominated — matter- 
dominated cosmological model. 

4. Theoretical predictions and experimental limits for stochastic gravita- 
tional waves. 

5. Variancies for coherent and squeezed quantum states. 

6. Classical trajectories at the compactified (j), <pj phase diagram. 

7. Classical paths in the (a, <fi) configuration space. 

8. The potential V(a) . 

9. Two branches of the caustic line and its singularity. 
10. Possible spectra of relic gravitons. 
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